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Abstract 

The main aim of the paper is to present the analytical solution of the Belavkin 
quantum filtering equation for damped harmonic oscillator being initially in the 
squeezed coherent state for diffusion observation with complex white noise. The 
comparison of the a priori and a posteriori mean value of the optical quadrature 
operators and the photon number operator is given. 

1 Introduction 

The quantum stochastic model of the Markovian dynamics of an open system under 
continuous nondemolition observation was given by Belavkin in pp. The solution of a 
quantum stochastic equation corresponding to an initial Gaussian state was described in 
detail by Belavkin in [1^ and for the general multidimensional case in ^ . The linear version 
of the Belavkin quantum filtering equation for the damped harmonic oscillator with the 
initial coherent state and a particular type of the Schrodinger cat state (superposition 
of two coherent states) was studied in more physical language in [3 E]. In this paper 
we give an analytical solution of the Belavkin quantum filtering equation for damped 
harmonic oscillator in the diffusive case with complex white noise and in the counting 
case for the system starting in a squeezed coherent state. We present the discussion on 
the eigenvalue equations for states preserved under stochastic time evolution. In the last 
section of the paper we investigate the expectation values of quadrature and the photon 
number operators of the system for a non-selective dynamics. 
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Poland, Electronic address: amdabrowska@gmail.com 
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2 The model 



We consider a simple quantum system S, namely a harmonic oscillator, interacting with 
the environment modelled by a Bose field given initially in the vacuum state. Instead of 
the usual description of a master equation for the Markovian approximation of this in- 
teraction we will start from a quantum stochastic Markovian model for unitary dynamics 
of the combined interacting system which gives a more detailed description for such phe- 
nomenological models and allows explicitly to consider not only the system but also the 
output fields introduced in pp. The unitary evolution operator for the system interacting 
with Bose field in the interaction picture with respect to the free dynamics of the field 
satisfies the quantum stochastic differential equation (QSDE) |2] 



dU{t) 



LL 1 

y^adB'^ (t) — ^/Jla) dB{t) — —a)adt — j^Hdt 



U{t), (1) 



with initial condition U (0) = /, where H = huo \ a)a + | j is the hamiltonian of the 
oscillator iS, a denotes the annihilation operator, is the adjoint creation operator and 
is a real coupling constant. The input field processes B^{t), B{t) describing the interaction 
up to the time t are related to creation and annihilation field operators &^(t), h{t) as follows 

t t 
B\t) = J b\s)ds, Bit) = J b{s)ds. (2) 



All information about S which can be realized in the course of the continuous nondemo- 
lition measurement up to time t is contained in the output field which is described by 
the Heinsenberg transformation of the input field: 

B{t) = U\t)B{t)U{t), B\t) = U\t)B\t)U{t). (3) 
Thus the photocounting signal can be expressed by the output counting process 

N{t) = U\t)N{t)U{t), (4) 

t 

where N{t) = J b^{s)b{s)ds. In the balanced homodyne detection when the output field 



is mixed with intense laser field at the same frequency (a local oscillator) we measure the 
observable |T2] 

dB\t)e''^ + dB{t)e-''^, (5) 

where is the phase of the local oscillator. When we assume that the system S is 
initially in a pure state denoted by \0{O)), the a posteriori unnormalized wave function 
\0{t)) satisfies the linear Belavkin filtering equation which for the counting process has 
the form [S] 

d\m = - - \^(i))dt \m)dNit). (6) 
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The equation is understood in the Ito sense, in which case 



K = ^H+^a^a. 
h 2 

When the system S is coupled to two independent modes of the Bose field we can get the 
diffusion observation with complex white noise. In this case the QSDE for the unitary 
evolution operator can be written in the form [H] 



dU{t) 



Uit), (7) 



where A{t) = (-Bi(t) + iB2{t)) /\/2 is the combination of two independent modes of the 



field and relevant coupling operators are defined by \J /i/2a, and — /i/2a. We can reahze 
the continuous observation of two independent components of the output field process 
being the Heisenberg transform of the input field process 

W{t) = + B\{t)) + I (B^it) + B\{t))] = W^{t) + iW2{t) . 

Note that W{t) = Wi{t) + iW2{t) is a complex white noise and for its real and imaginary 
part one can write {dWi{t)f = {dW2it)f = \dt, dWi{t)dW2{t) = 0. 
The unnormalized wave function conditioned by the trajectory up to t of the observed 
output field process W{t) = W{t)W{t)U{t) satisfies the following equation 

d\^it)) = -[^H + |ata) \m)dt + ^a\m)dWit) . (8) 

When the results of measurement are not read, i.e. there is no selection of the trajectory, 
we get a priori dynamics which is described by the Markovian master equation of the 
Lindblad-Kossakowski form [3 |H] 



apa^ — - (a^ap + pa^aj . (9) 



3 Analytical solutions for the diffusion observation 
with complex white noise 

Let us first describe the solution of Eq. © for |(^(0)) = \a), where \a) is a coherent state 
and a G C is the complex amplitude of this state. In this case the normalized posterior 
state \ipit)) = II \(pit)) II coincides up to a phase with the prior coherent state \a{t)) 

as proved by Barchielli and Belavkin p for the more general case of a single-mode field 
in a cavity with some source interacting with thermal bath. Since the Eq. (jH)) preserves 
the coherent state we can insert \(p{t)) = g{t)\a{t)). This provides us with two differential 
equations: 

da{t) = -(^iuj + a{t)dt, (10) 
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dg{t) 



-^-^dt + ^a{t)dW{t) + ^d (la(t)l' 



9{t), 



with initial conditions g{0) = 1, a{0) = a, and hence 



\0{t)) = exp 



t 



(12) 



where a{t) = ae"(^"+^)*. 

The expectation value of any system operator Z conditioned by a trajectory of the results 
of continuous measurement is given by 



(13) 



imim) 

This formula defines the a posteriori mean values of quadrature operators X = ^ (^a + j , 

{X{t)) = Rea{t), {Y{t)) = Ima(t), (14) 
and of the photon number operator n = a^a 



{h{t)) = Ht)f 



a\ 



(15) 



One can notice that the obtained values do not depend on the measured noise. The 
linearity of Eq. (jH)) allows us to find easily an exact solution for any superposition of 
coherent states. Let us focus on the results for two initial states of the form 



l^(0)) = |«,±) 



|q;) ± I — a) 



N. 



± 



(16) 



where N± = y2(^l±e jg normalization. The state |a, +) is called the even 

coherent state because it is a superposition of states with even numbers of photons, 
and accordingly, the state \a, — ) is called the odd coherent state. When \a\ ^ 1 the 
components | ± a) become macroscopically distinguishable and we get two examples of 
the so-called Schrodinger cat state. The analytical solutions for (fTB|) are given by the 
expression 



\0{t)) = K{t) (e^W|aW) ± e"^^*^| - « 



;i7) 



iujt I 1 



+ i|ar(e-^*-l) /N± and a{t) 



where x(^) = / a{s)dW{s), K{t) = exp 



^g-(ja;+^)t^ We can assume that for large values of amplitudes |a| 1 and time t ^ l//i 



{a{t)\-a{t)) 



0. 



;i8) 



Hence, for the even coherent state one can obtain the formulae 
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(X(t)) = Rea{t) tanh (2Rex(t)) , {Y{t)) = lma{t) tanh (2Rex(t)) (19) 
and {n{t)) = e^^* |«|^. The uncertainties of the quadrature components have the form 

AX{t) = J- + Re^a(t) (l - tanh^ (2Rex(t))), (20) 



AY{t) = y ^ + lm^a{t) (l - tanh^ (2Rex(t))). (21) 

The expression tanh (2Rex(t)) tends to assume one of two values ±1 which means that 
one can observe the selection one of two solutions for two considered coherent states. We 
note that for the admitted assumption only {n(t)) does not depend on the noise terms. 
For others details of the solution of the quantum filtering equation for the Schrodinger 
cat state see jl]. 

Now we investigate the other analytical solution of Eq. (jSJ for the case when the system 
is initially in a squeezed coherent state defined by the formula 

1^(0)) = S{OD{a)\0) = S{0\a) = (22) 

where 

^(0 = expQeV-^^(at)'), (23) 
is known as the squeeze operator (with ^ = ge''^ G C), and 

D{a) = e"^^-"*^, (24) 

is called the displacement operator. There exists an equivalent definition of a squeezed 
coherent state \p, = D (p) S (0 |0) with D (/3) S {^) = S (0 D (a), where [HI 

a = P cosh g + /?*e*^ sinh g, (3 = a cosh g — a*e^^ sinh g. (25) 

We can use the above relations to transform the mean values of the operators for one of 
these two states to the mean values for the second one. 

The solution corresponding to chosen initial state (j22p can be written in the form 



\m) = G{t)sm)Ht)), 



(26) 



which simply means that squeezed coherent state is preserved under stochastic evolution. 
The proof of it was given by Belavkin in [U^ . We insert (j^UI) in Eq. (jH} and multiply the 
equation from the left hand side by S'^{C,{t)). To find the increment dS {^{t)) we use the 
form of the squeeze operator with normal order of annihilation and creation operators 

uni 



S{0 = (cosh g) 



-1/2 



exp 



In (cosh g) a' a 



e 2 



(27) 
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where F = e*^ tanh q. The Baker-Hausdorff formula allows us to obtain unitary transfor- 
mations 



S\^)aS{S,) = a cosh q — a)e^^ sinh q, 
S\i)a)S{i) = cosh g - ae''^ sinh 



(28) 
(29) 



S^{i) {a))^ S{i) = (at)%osh2^- (20^0 + l)e-^^sinh^ cosh ^ + aV2^^sinh2^, (30) 



S\i)a)aS{i) = a) a cosh^ q + (a^a + l) sinh^ q - a^e"^^ + (a)) 
These formulae allow us to obtain the differential equations: 



sinh ^ cosh ^. (31) 





d9{t) 


= -2ujdt, dg{t) = - 


-/I sinh ^(i) cosh g{t)dt , 


(32) 


da{t) = 


— (^iuj 


+ -/xsinh^^(t) 


a(t)dt - ^e'^(*) sinh g(t)dW(t) , 


(33) 


dG{t) 
Git) 


icu , 

dt- 

2 


f ^d \a{t) f + ^/Jla{t) cosh g{t)dW{t) - ^ sinh^ g{t)dt+ 





+ iie~'^^^'^a'^{t) sinh ^(t) cosh g{t)dt , 



(34) 



with the initial conditions G(0) = 1, a{0) = a, 9{0) — 9 and ^(0) = g. Integrating these 
equations yields 



9{t) = 9- 2ujt, g{t) = arc tanh (e"''* tanh g) , 
a{t) = e-(-+t)* [a - ^e^' smh g } e-i-^^>dW{s)\ , 



(35) 
(36) 



^ /COS11^^_^+1(|,(,)|2_|^(0),. 

V cosh^ 



X 



X exp 



J (a{s) cosh g{s)dW{s) + -^JJie '^^*^q;^(s) sinh g{s) cosh g{s)ds^ 



Now one can find the a posteriori mean values of quadrature operators 
{X(t)) = Re (a (t) cosh ^(i) - a* (t) e^^^*) sinh ^(i)) , 



(37) 



(38) 



{Y{t)) = Im cosh g{t) + a(t)e^^^W sinh g{t)^ . 



(39) 



In contract to the analogous formulae for the coherent states, these expressions depend 
on the noise terms. The uncertainties of quadrature operators for considered squeezed 
coherent state have the form 



AA'(i) = - 



Am = I 



1 + 



1 + 



2e tanh g 
1 - e-^A** tanh^ g 

2e~^* tanh g 



e tanh g — cos I 



e ^* tanh g + 



g + cos I 



2uot) 



+ 2ut) 



1/2 



nl/2 



1 - e-2Mt tanh^ g 

and one can check that for fi = (unobserved system) we get expressions 



AX{t) 



1 



1 + 2 sinh g [sinh g — cosh g cos {6 — 2ujt)], 



(40) 
(41) 

(42) 



AY{t) = ^ y'l + 2 sinh g [sinh g + cosh g cos (61 - 2cji(;)] . 



(43) 



Thus the squeezing determined by coefficient g{t) decreases in time. We get the squeezing 
in the X-quadrature whenever cos [6 + 2ut) > e~^* tanh g and in the ^-quadrature when 
cos {9 + 2ojt) < -e"^* tanh g. 

One can note that when we put = in the above formulae we get the solution for the 
initial coherent state and when the amplitude a = then all expressions correspond to 
the initial squeezed vacuum state defined by S'(^)|0). 

For all three mentioned initial states the uncertainties of quadrature operators are inde- 
pendent of the noise terms. Obviously any a posteriori and a priori state tends asymp- 
totically to the vacuum because we consider an undriven open harmonic oscillator. 



4 Conditional states in the counting process 

In this section we will discuss direct detection of the output field 13 {t) and B\t) given by 
The evolution between the counts (i.e. when dN{t) = in Eq. (jU))) can be written 
for unnormalized wave function in the form 

^^ = -AW)). (44) 

When a count of a photon is registered at time t = to {dN(to) = 1 in Eq. © and the 
other term in the right hand side of Eq. (jH)) is neglected) the normalized state of the 
system immediately after the count reads 

|#o + dt))^ , "l-^"'" . (45) 

{(p{to)\a^a\(p{to)) 
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We should recall that for simplicity we abandon the retardation due to the time of flight 
of the photons to detector. The solution of Eq. (|44|) can be written in the form 



\0{t)) = exp 



iut 

— - — + 



a^at 



(46) 



where |<^(0)) denotes initial state of the system. Therefore for the initial coherent state 
I a) we get 



\0{t)) = exp 



iut \a\ / 
9. 9~ V ~ 



-fit 



\a{t)), 



(47) 



and hence (X(t)) = Rea(t), {Y{t)) = Ima(t) and (n(t)) 



-fit 



a\ 



where a{t) 



The probability of no counts up to time t is given by P{t) = II |^(t)) IP 0. For the 
solution (P7j) we obtain 



P{t) = exp 



\a\ 



(1 



-fit 



(4J 



which in the limit t ^ oo gives the expected probability to be in the vacuum for the initial 
coherent state. One can also easily find a similar formula for the initial even coherent 
state 

(49) 



P{t) = cosh (lal^e-'^*) / cosh 



a\ 



and lim P{t) = 1/ cosh |a|^. For the initial odd coherent state 



t— >oo 



P{t) = sinh(|a|^e"''*) /sinh 



a\ 



(50) 



and lim P{t) = which implies that for the initial odd coherent state there is at least 

t^oc 

one count in the time interval (0,+oo). The a posteriori mean values of quadrature 
operators for both initial states are equal zero and the expectation values of number 
operator changes in time according to the formula 



{nit)) 



|a|^e-'^*tanh(|a|^e-'^* 



for the even coherent state, while for the odd coherent state 



{n{t)) 



coth (^1 



(51) 



(52) 



In the limit t ^ oo, they give us expected value zero. For the uncertainties of quadrature 
operators we get the damped oscillating functions of the form 



AX(t) 



1 + ARe^ait) + (l - 41m^a(t)) e^^l" 



|2g-Mt 



1/2 



2Vl + e-2|"l 



(53) 
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AY{t) - 
for |<^(0)) = |a;,+), and 

AX{t) -- 
AYit) 



1 + 4Im^a{t) + ( 1 - 4Re^a(t)) e 



1/2 



2V1 + e-2|"re- 



(54) 



1 + ARe^a{t) - (l - 4Im2a(t)) e'^l"!'--''* 



1/2 



2\/l 



-2|ape-M* 



1 + ilm^ait) - (1 - 4Re^a(t) ) e 



1/2 



2\/l 



(55) 



(56) 



for |<^(0)) = ). One can easily check that for the considered case the state of the 
system after a count at time t = to becomes 



IHto + dt)) 



e 2 



a 


a{to)) T 




a 


(to)) 






-2 


a 


2 Q-^J,to^ 



(57) 



which means that the harmonic oscillator jumps from one type of "cat" to another type 
of "cat" 0. 

Now let us now briefly present the fate of squeezed coherent state when there are no 
counts. To find the solution of Eq. pij) for |<^(0)) = \^,a) we can use the insertion (|^. 
This procedure provides us with solutions 



cosh^(t) _/ 

a[t) = fi V 

cosh g 



(5^ 



G{t) 



X exp 



V cosh g 
(^cosh^ g — sinh^ ge~'^^^ — 1^ 



ia^g-i^'W coth g 5 5 ^ 

2 cosh g — sinh ge ^m* 



(59) 



where G(0) = 1, a(0) = a, 6'(0) = 6*, ^)(0) = g and the formulae for 6 (t) and (t) coincide 
with expressions (jH^j) . We can check that P(t) = |G'(t)|^ gives us the known limit ^1] 



lim Pit) = e-l°l' (cosh g)~^ exp fee f e^^a^) tanh p 

t— >oo L V / 



(60) 



the probability of being in the vacuum for the initial squeezed coherent state |,^, a). 
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5 The eigenvalue equations for the preserved states 

We show a technically easier way of obtaining differential equations for these coefficients 
which allow us to get the a posteriori mean values of the operators. Since a coherent state 
is preserved under stochastic evolution it follows that if the system is in a coherent state 
at instant t, 

a\m) = (^it)\m), (61) 

at instant t + dt we get 

a\(f{t + dt)) = a{t + dt)\0{t + dt)). (62) 
From these equations one obtains 

(a - a{t) - da{t)) d\ip{t)) = da{t)\(p{t)) , (63) 

where d\if{t)) = \0{t + dt)) — \0(t)). Inserting d\(f{t)) from the Eq. © or Eq. (jHI) into 

Eq. ()63|). we can obtain the relevant differential equation for a (t). 

The eigenvalue equation for the considered squeezed coherent state has the form 

S{OaS\0\^,a) = a\^,a). (64) 
Rewriting this equation in a more convenient form by putting 

S{OaS\0 = + a^T2 , (65) 

where Fi = cosh g, F2 = e*^ sinh g, and applying similar reasoning to that for the coherent 
state we get 

a (Fi(t) + dFi(t)) + at (F2(t) + rfF2(t)) - a(t) - rfa(t)] d\^{t)) = 

= (c/a(t) - adTi{t) - aUV^if)) \(p{t)) . (66) 

This provides us with the relevant differential equations for the coefficients Q if)-, g it) and 
a{t). 



6 The comparison of the a priori and a posteriori mean 
values operators 

We now investigate the a priori mean values of the operators. The nonselective dynamics 
given by the dynamical semigroup takes in the considered case the form Eq. Q. We 
do not need to solve this equation for studying the required characteristics. The a priori 
mean value of the number operator is defined as 

{n{t)) = Tr (atap(t)) , (67) 
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hence obviously 

^ = Tr(atap(t)), (68) 
One can check that for p(t) given by Eq. Q the equation (|UHj) leads to a simple formula 

{n{t)) = e-^'{n{0)). (69) 
In analogous manner we can derive 

{X{t)) = [{X{0)) cosujt + (r(0)) sinujt] , (70) 

{Y{t)) = e'^ [{Y{<d)) cosujt - (X(0)) sinut] . (71) 

Thus one can notice that for the coherent state the a priori and a posteriori mean values 
of operators of the system coincide. If the system is initially in the even coherent state 
we get for the photon number operator the mean value of the form 

(n(t)) = e"'^* |a|^tanh|a|% (72) 
and for the odd coherent initial state 

{n{t)) = e-^'\afcoth\a\\ (73) 

Therefore for |a| ^ 1 and t -C l//i in the diffusive case the conditioned and unconditioned 
mean values of the operator n are undistinguished in contradiction to the expectation 
values of quadrature operators which in semigroup dynamics do not depend on time and 
remain equal zero. One can also check that the squeezed vacuum state has only nonzero 
expectation value of the photon number operator 

(n(t)) = e-'^^sinh^^. (74) 

Because (X(0)) = (X(0)) and {Y{0)) = {Y{0)), we can easily find the a priori mean 
values of the operators for the squeezed coherent state using the expressions (jHHj) and 
(jSni)- To obtain formula for the mean value of the photon number operator it is sufficient 
to know that for this state 

(n(0)) = \a cosh g — a* sinh gf + sinh^ g. (75) 

7 Conclusions 

We have obtained the analytical solution of the Belavkin quantum filtering equation for 
damped harmonic oscillator starting in squeezed coherent state making use of the fact that 
the considered state is preserved under stochastic evolution. From the solution for the 
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squeezed coherent state we can get the expressions describing the fate of initial coherent 

state and squeezed vacuum state. 

It was noted that for the undriven oscillator being initialy in coherent state the a priori 
and a posteriori states coincide up to the phase. We have checked that the uncertainties 
of quadrature operators for squeezed coherent state in diffusion observation have the form 
of damped oscillating functions independent of the noise terms and that the parameter 
g{t) describing the squeezing changes in time in a smooth non-random way. 

8 Acknowledgment 

1 am grateful to Prof. V. P. Belavkin for his hospitality and discussions, and 1 am indebted 
to Prof. P. Staszewski for careful reading of the manuscript. This work was supported by 
the EU RTN grant "QP and Applications" . 

References 

[1] V. p. Belavkin, Optimal Measurement and Control in Quantum Dynamical Systems, 
Preprint No. 411, Institute of Physics, Nicholas Copernicus University, Toruh, 1979. 



[2] 


V. 


P. Belavkin, Rep. Math. Phys. 43, 405 (1999). 




[3] 


p. 


Goetsch and R. Graham, Phys. Rev. A 50, 5242 (1994). 




[4] 


p. 


Goetsch, R. Graham and F. Haake, Phys. Rev. A 51 136 (1995) 




[5] 


A. 


BarchieUi and V.P. Belavkin, J. Phys. A 24, 1495 (1991). 




[6] 


V. 


P. Belavkin, Commun. Math. Phys. 146, 611 (1992). 




[7] 


G. 


Lindblad, Commun. Math. Phys. 48, 119 (1976). 




[8] 


V. 


Gorini, A. Kossakowski and E. C. G. Sudarshan, J. Math. Phys. 


17, 821 (1976). 


[9] 


B. 


M. Garraway and P. L. Knight, Phys. Rev. A 50, 2548 (1994). 




[10] 


See for example: D. F. Walls and G. J. Milburn, Quantum Optics 


(Springer Study 



Edition, Heidelberg, 1994). 

[11] See for example: L. Mandel and E. Wolf, Quantum Coherence and Quantum Optics 
(Cambridge University Press, Cambridge, 1995). 

[12] See for example: M. O. Scully and M.S. Zubairy, Quantum optics (Cambridge Uni- 
versity Press, Cambridge, 2001). 



12 



